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ABSTRACT 


Studies  on  optimization  of  a  single-stage  rocket  perform¬ 
ance  with  and  without  the  use  of  the  Vacuum-Air  Missile  Boost 
(VAMB)  system  are  presented  here.  Optimization  is  defined  as 
maximizing  the  altitude  to  which  a  certain  rocket  of  a  given  weight 
and  propellant  can  send  a  fixed  payload.  Standard  atmosphere  has 
been  used  for  the  analysis.  Performance  estimation  and  optimiza¬ 
tion  of  a  hypothetical  1,000  lb  rocket  with  and  without  the  VAMB 
system  have  been  done  under  the  simplifying  assumptions  of  con¬ 
stant  (neglecting  the  pressure  thrust)  for  the  powered  flight  and 
a  constant  drag  coefficient  for  the  whole  flight.  Similar  studies 
have  also  been  performed  for  eighteen  real  rockets  without  the 
VAMB  system.  It  was  found  that  the  optimum  thrust-to-weight  ratio 
(T/W)  for  conventional  rocket  can  be  estimated  roughly  for  prelim¬ 
inary  design  of  rockets  under  these  assumptions.  The  maximum 
altitude  reached  is,  however,  underestimated  by  this  approach. 

For  launch  of  rockets  utilizing  the  VAMB  system,  determination 
of  optimum  T/W  is  also  important.  It  becomes  critical  when  the 
loading  density  of  the  rocket  is  low  (£.10)  ,  its  fuel  ratio  ap¬ 
proaches  unity  and  the  initial  velocity  imparted  is  high  (£500  fps). 
Detailed  performance  and  optimization  studies  have  been  presented 


on  the  Areas  rocket.  These  are  based  on  measured  data  of  thrust 
and  drag  coefficient  variations  for  a  recent  Areas  design,  as  ob¬ 
tained  from  the  White  Sands  Missile  Range.  It  was  found  that  the 
optimum  T/W  for  Areas  with  and  without  VAMB  system  is  approxi¬ 
mately  2.4,  which  is  almost  half  of  the  present  design  value  {4.5). 
The  gain  in  maximum  altitude  of  the  optimized  Areas  launched  with 
an  initial  velocity  of  1,000  fps,  over  that  of  the  present  design 
Areas  without  booster,  was  estimated  to  be  nearly  93  per  cent. 

It  is  concluded  that  the  VAMB  system  is  a  very  helpfui  adjunct 
to  small  rockets  of  optimized  design. 
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AND  OTHER  SINGLE-STAGE  ROCKETS  KITH 
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The  vacuum-sir  missile  boost  systems  was  developed 
earlier :  * z  * 3  * **  * £  * 7  for  the  boosting  of  small  to  medium  size 
solid  propellant  rockets,  by  imparting  them  an  initial  velocity 
at  launch.  Host  of  the  previous  work  was  concerned  with  the  flow 
mechanics  and  dynamics  of  the  missile  inside  the  tube.  -Thus, 
reliable  estimates  of  the  velocities,  which  can  be  imparted  to 
a  certain  rocket  by  this  system,  can  be  made.  The  problem,  of 
optimising  the  rocket  performance  by  utilising ' this  system, 
however,  remains  unsolved.  This  thesis  essentially  studies 
some  of  these  aspects  with  special  reference  to  a  particular 
rocket,  the.  Areas  (All  Purpose  Rocket  for  Collecting  Atmospheric 
Soundings) _ 

The  vacuum-air  missile  boost  system,  later  referred  a 
as  the-  VAMB  system,  consists  of  a  partially  evacuated  vertical 
launching  tube  with  a  breakable  seal  at  the  top  and  with  the 
lower  end  sealed  by  the  missile  on  a  held-down  sabot.  When  the 
sabot  is  releasee,  atmospheric  pressure  from  air, entering  the 
base  of  the  tube  accelerates  the  missile.  A  muzzle  chamber  near 
the  too  of  the  tube  reduces  the  velocity  loss  due  to  compression 
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of  the  residual  air  above  it,  and  the  tube  length  bey-nd  the 
muzzle  chamber  serves  as  an  atmospheric  shock  reducer.  The 
velocity  imparted  by  this  system  can  be  up  to  the  sonic  veloc¬ 
ity.  The  VAMB  system  is  considered  suitable  for  upper-atmcs- 
pheric  sounding  rockets.  Some  of  the  advantages  of  the  system 
when  compared  to  the  usual  surface-launching  technique  are  a 
significant  fuel  saving  or  an  increase  in  payload5,  less  sensi¬ 
tivity  to  surface  gusts,  and  a  shelter  prior  to “firing. 

Much  work  has  been  done  oh  the  optimum  performance  of 
the  usual  surface-launching  rockets’ • 1  0 • 1 x» 1 2 .  For  optimization 
of  the  performance  of  a  rocket  launched  by  the  VAMB  system,  the 
present  study  was  initiated  to  investigate  the  effects  of  the 
various  rocket  parameters .  Variations  of  parameters ,  such  as 
initial  thrust- to-weight  ratio  and  fuel  ratio,  on  the  perform¬ 
ance  of  sounding  rockets  have  been  studied.  Essentially,  the 
main  problem  considered  herein  ds'  that  of  determining  rocket 
burning  programs  that  use  a  given  amount  of  fuel  to  propel  a 
given  mass  to  a  maximum  altitude.  Or  putting  it.  differently, 
the  criterion  of  optimization,  for  this,  investigation,  is  to 
maximize  the  altitude  to  which  a  given  payload  can  be  sent  by 
a  rocket  of  a  specified  weight. 

Preliminary  study  was  conducted  with  a  hypothetical 
rocket  of  1,000  lb.  with  characteristics  similar  to  comparable 
weight  rockets,  JL  wide  variety  of  real  rockets  were  then 
investigated.  Further  refined  study  was  made  specially  for 
the  single-stage  Areas  rocket  currently  in  use.  Finally, 
forward  stagnation  temperatures  and  .pressures  of  the  Areas 
rocket  during  optimum  performance  were  evaluated. 
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Its  of  this  investigation  indicate  that  a 

-gain  in  maximum  altitude  can  be  attained  by  simul- 

mization  of  the  thrust-to-weight.  ratio  and  applica- 


CHAPTER  II 


FORMULA? ION  OF  THE  EQUATIONS  OF  MOTION 


2.1  Simplified  Equations  of  Motion  for 
Vertical  Trajectory  . 


In  developing  the,  equations  of  motion  for  the  vertical 
flight  of  a  single-stage  rocket  in  an  atmosphere  near  the 
earth’s  surface,  the  following  assumptions  are  made: 

i)  The  earth's  gravitational  field  is  constant  and 
the  gravitational  acceleration  is  taken  to  be  32.2  ft/sec*,  the 
val>  at  the  sea-level.  Although  it  decreases  with:  increasing 
altitude,  the  error  induced  is  about  10  per  cent -at  the  alti¬ 
tude  of  200  miles  which  is  the  extreme  altitude  within  our  scope 
of  study. 

ii)  The  density  and  pressure -of  the  atmosphere  vary 
parabolically  with  X he  altitude  up  to  55,000  ft.  and  expor.en-  - 
tially  above  this  altitude.  The  expressions  fc-r  the  variations 
fit  the  experimental  data  fairly  closely  as  described  in  Sec*  3.2. 

iii)  The”  rocket  propellant  is  consumed  at  a  constant 
rate.  This  is  generally  accepted  by  rocket  designers. 

iv)  The  exhaust  velocity  of  the  burned  gases  relative 
to  the  rocket  is  constant-  This  is  quite  reasonable  for  most 
rockets . 

(5) 
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v}  The  thrust  direction  is  aligned  with  the  rocket 
longitudinal  axis  which  a 1 wavs  remains  vertical.  For  vertical 
flight,  this  assumption  is  quite  acceptable. 


the  earth's  rotation  ere  neglected.  This 


is  to  be  justified  in  tae  following  section. 

viz)  Effects  of  the  atmospheric  winds  are  neglected. 

The  wind  speed  is  difficult  to  estimate  and  its  magnitude  is 
much  smaller  than  the  velocity  of  the  rocket  thus  having  little 
effect  on  the  trajectory. 

vlii)  Only  thrust,  aerodynamic  drag  force,  and.  gravity  - 
affect  the  motion  since  other  forces  are  small  compared  to  these 
ix)  The  aerodynamic  drag  force  is  given  by 


P  =  7  CDpAV2 


where  A,  V  are  the  trontal  area  and  the  velocity  of  the  rocket; 
Cr>  is  the  drag,  coefficient -which  may  be  a.  constant  or  a  variable 
{a  function  of  Mach  number);  and"  p,  the  ambient  air  density. in 
slug/ ft 3 ,  is  a  function  of  altitude. 

Under  these  assumptions  the  equation  of  motion  for  a 
vertical  trajectory  during  burning  (sustain  phase)  is  written  as 


{2-2} 


T  =  thrust  =  W  Vg/g  *  Ae  {PS£,  -  P£} 

D  =  aerodynamic  drag  =  *  Cq  pAV2 
W  =  lift-off  weight 
v:  -  propellant  flow  rate  in  Ib^/sec 
t  -  time 

"g  =  gravitational  acceleration 
Vg  =  exhaust  velocity.  of  ejected  gas 
Aq  =  effective  area  of  nozzle  exit 
PsX.  =  atmospheric  pressure  at  sea-level 
Pa  =  ambient  pressure  of  atmosphere 
Also  by  definition 


where 

Y  -raititude-  in"  feet. 


(2-3) 


-  (2-4) 


In  Sg.  {2-2} ,  the  thrust  includes  pressure  thrust, 

Ae  {Pgj^  -  Pa)  ,  in  addition  to  momentum  thrust,  flWg/g.  The 
pressure- -thrust,  usually  not  so  important,  is  only  a  few  per 
cent  of  the  momentum  thrust  and'  hence ,  is  neglected  in  the 
preliminary  study. 

Once  burning  is  terminated,  only  aerodynamic  drag  and 
gravity  ^affect  the  free  flight  (coasting  phase)  of  the  rocket 
and  -the  equation  of  motion  reduces  to 


where  Wp  is  the  weight  of  the  propellant. 


In  general,  a  x  a  x  £  is  small  compared  to  u  x  ts  x  R, 
and  thus  may  be  neglected.  This  is  based  on  the  argument  that- 
S  is  of  the  order  of  200  miles  or  less  while  R  is  of  the  order 
of  4,000  miles.  Furthermore,  one  may,  with  little  loss  of 
accuracy,  assume  that  the  external  forces  -  thrust,  drag  and 
gravity  -  are  acting  in  the  direction  of  local  vertical,  i.e., 
along  the  y-axis.  The  equation  of  motion  in  the  inertial  coor¬ 
dinate  system  is  then  written  as 

(T  -  D  -  mg)  J  =  mj^a  x  a  x  R  +  2a  x  S 
Introducing  a  =  acos  4>  X  +  a  sin  6  J, 

one  obtains  the  scalar  form  of  the  equations  of  motion  in  the 
coordinate  system  xyz  as  follows: 


x  : 

0  =  mja2  R  sin  $  cos  6+2 

9 

az 

sin  6  +  jij 

(2-11) 

y  : 

T  -  D-mg  =  m  ^-a2  R  cos2$  - 

-  2 

az  cos  6  +  Y  ] 

(2-12) 

z  : 

0  =  ffi  £ 2u  (y  cos  $  -  x  sin 

*> 

+  z  1 

(2-13) 

+  S 


xyz 


(2-10) 


xyz 


The  rocket  is  travelling  nearly  vertically  upwards, 
so  x  «  y  and  i  <<  y,  and  the  products  ax  and  a£  are  vanishingl 
small  because  of  the  smallness  of  u>{0.727  x  10-1’  rad/sec).  Thu 
the  Eqs.  (2-11)  through  (2-13)  reduce  to 
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vJ  *”  Iv  X  ^ 


:os  q  +  x 


(2-14) 


_  pi  _  _  = 


CO$  S  -5- 


(2-15) 


=  ,  9  \T  - O 


(2-16) 


2n  examining  Sc*  (2-155  /  the.  centrifugal  acceleration, 

R  -  0.11  ft/sec£,  is  only  G.CG3g  and  thus  may  well  be  neglected 
Since,  even  for  a  long  time  of  flight,  say  200  seconds,  its 
contribution  to  the  altitude  only  reaches  2,200  feet  which  i3 
no  more  than  one  per  cent  of  the  maximum  altitude.  Therefore, 

£q.  (2-15)  reduces  to  the  same  form  as  derived  by  assuming  non¬ 
rotating  earth.  This  justifies  the  assumption  that  the  earth's 
rotation  has  little  effect,  on  the  maximum  altitude  of  the  mis¬ 
sile  within  our  scope  of  study.  : 

It  can  be  seen  from  Has.  (2-14)  and  (2-15)  that  the 
deflections  in  the  x  and  z  directions  depend  on  the  latitude 
of  the  launch  spot,  4. 

In  the  x-cirection,  there  would  be  no  deflection  if 
launching  is  done  at  the  equator  or  at  the  poles.  If  g  =  45°, 
however,  it  reaches  the  maximum  and  Eq.  (2-14)  becomes 


~v , 


Upon  integration  one  obtains 


x  =  -0.25 
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x  =  -0 . Q276t2  {2-17} 

For  various  times  of  flight,  the  deflections  in  the  x-direction 
for  $  =  45°  are  indicated  below  {Table  2-1)  . 

The  deflection  in  the  z-direction  vanishes  as  j  =  S0°r 
while  it  reaches  its  maximum  as  C  =  0°  (at  the  equator) .  It 
can  be  estimated  by  taking  an  average  velocity,  yav, during  the 
flight  so  that  Eq.  (2-16)  may  readily  be  integrated.  By  doing 
so,  one  obtains 

z  =  -  ^yav  cose  t2  (2-18) 

The  deflections  in  the  z-direction  for  various  times  cf 
flight  and  average  velocities  are  presented  in  Table  2-2. 

The  simplified  analysis  of  the  earth's  rotational  effects, 
which  has  so  far  been  conducted,  indicates  that  the  earth's 
rotation  has  little  effect  on  the  maximum  altitude  of  the  missile 
at  least  within  our  scope  of  study,  and  hence  has  no  significant 
contribution  to  The  t  ^flections  of  the  missile  in  the  x 

and  z  directions,  however,  can  be  of  considerable  significance 
for  long  time  cf  flight  in  certain  launching  areas  as  indicated 


in  Tables  2-1  and  2-2. 
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In  determining  the  optimum  thrust- to -weight  ratio  for  a 
specified  rocket  with  an  initial  velocity  imparted  by  the  VAM5 
system,  rocket  parameters  such  as  lift-off  weight,  specific 
impulse,  fuel  ratio  (ratio  of  weight  of  propellant  to  lift-off 
weight),  and  frontal  area  are  fixed  while  the  thrust- to-weight 
ratio  is  varied.  Change  c-f  thrust- to-weight  ratio  is  related 
to  a  definite  change  in  burning  time  or  flow  rate,  area  of  nozzle 
exit,  and  curve  for  variations-  of  drag  coefficient  with  Mach 
number  durii^g  burn  phase . 

Kith  some  specified  thrust- to-weight  ratio  and  rocket 
characteristics,  one  can  readily  compute  the  weight  of  propel¬ 
lant,  propellant  flew  rate,  and  burning  time.  Having  this 
information,  the  equations  of  motion  can  be  integrated  numeri¬ 
cally  to  obtain  a  value  for  maximum  altitude.  r 

One  then  changes  the  thrust-to-veight  ratio;  and  computes 
the  new  values  of  flow  rate,  burning  time,  area  of  nozzle  exit- 
and  Cq  curve  during  burn  phase.  The  equations  of  notion  are 
integrated  again  to  produce  a  new  value  of  maximum  altitude. 

This  procedure  is  repeated  for  a  selected  range  of  thrus*:- to- 
weight  ratios.  A  curve  is  thus  plotted  for  maximum  altitude 
versus  thrust- to-weight  ra*-ic  and  the  peak  on  the  curve  locatec 
the  optimum  thrust-to-veight  ratio  for  the  specific  rocket  under 
the  specified  launching  conditions. 

3efore  numerical  analyses  proceed,  the  definition  of 
thrust-to-veight  ratio  should  be  identified.  It  is  defined  as 
the  ratio  of  average  sea-level  thrust  (assuming  optimum  .iozzle 
for  sea-level  operation)  to  lift-off  weight.  Under  this  defi¬ 
nition,  the  thrust  only  includes  the  momentum  thrust.  Thus, 


IS 

for  altitudes  of  55,000  ft.  or  less,  and  an  exponential  form 


Pi  e 


aY 


(3-2) 


for  altitudes  above  55,000  ft. 

.  The  constants  obtained  from  a  least  squares  curve  fit 
program  are 

Cc  =  0.757878  x  1G**1 
'  C,  =  -  0.203365  x  10"s 

C2  =  0.150855  x  10"10 
Cj  =  0.125184 

a  -  -  0.480263  x  10”1,  '  , 

These  constants  will  give  the  atmospheric  density  as  a  function 
of  altitude  in  lbm/ft3 .  _ 

Similar  approach  was  applied  to  variations  of  the 
atmospheric  pressure  in  lb/ft2  and  the  resulting  expressions 
are  _  _ 

P  s  ao  +  as  V  +  aj  {3-3} 

for  Y  <  55,000  ft,  and 

P  =  Pi  eYY  {3-4) 


for  Y  >  55,000  ft. 


The  constants  were  found  as  follows: 
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P  RE 


For  220,000  ft  <  Y  <  164,000  ft, 

T.  =  450  +  .JJQ —  (Y  -  120,000)  (3-6) 

a  44,000 

For  164,000  ft  <  Y  <  200,000  ft, 

T  =  630  °R  (3-7) 

cl 

And,  for  200,000  ft  <  Y  <  250,000  ft, 

T  =  450  +  -iPP- . .  (250,000  -  Y)  (3-8) 

a  50,000 

Since  experimental  data  for  temperature  are  not  available  for 
altitudes  higher  than  250,000  ft,  it  was  considered  necessary 
to  resort  to  the  perfect  gas  law  to  estimate  the  temperature 
above  this  altitude. 


:l  of  various  parameters  cn  the 


PERFORMANCE  01 
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KZ-.-iA  KwV, 


KETS  UNDE; 


4.1  Introductory  Comments 

As  mentioned  before,  the  complexity  of  the  equations  of 
motion  makes  it  impossible  to  obtain  an  analytical  mathematical 
expression  for  the  maximum  altitude.  Accurate  solutions  are 
obtained  only  by  using  a  digital  computer,  but  each  of  these 
applies  to  a  certain  rocket  under  some  specified  launching  con¬ 
ditions.  However,  by  examining  performances  of  many  rockets, 
the  tentative  influence  of  rocket  parameters  still  can  be  re¬ 
vealed.  It  is  the  purpose  of  this  chapter  to  investigate  a 
variety  of  real  rockets  to  acquire  some  general  ideas  concerninc 


:ne  performance  or  tne  rocket. 

For  simplifying  the  analysis 
lh  this  oreliminarv  study; 


jsumotions  are  maae 


5  v'jr?  i  f* 


ls  suited  to  be  constant  throuc 


ISC.  uU6  t-O  Cr.G  V  ci  IT  a.  d  -lOH  O I  Si  U I  t.Uut 


Tne  pressure 
r.eqlected  since  it 


.s  on iv  a  re*  oe: 


ii)  An  estimated  constant  drag  coefficient  is  assumed 


throughout  the  vertical  ascent  of  the  rocket. 

Although  these  assumptions  are  not  valid  for  the  actual 
performance  of  the  rocket,  results  are  yet  considered  helpful 
in  parametric  studies. 

Numerical  integration  with  time  step  of  0.5  seconds 
has  been  conducted  with  the  aid  of  a  digital  computer.  First 
of  all,  a  hypothetical  rocket  of  1,000  lb  is  investigated,  ar.d 
then  a  number  of  real  rockets  are  studied. 


4.2  Hypothetical  rocket  of  1,00'  lb 

The  Aerobee  rocket  was  the  first  rocket  of  its  type 
developed  specifically  to  investigate  the  upper  atmosphere  ar.d 
was  for  a  long  time  practically  the  only  readily  available 
upper-atmospheric  sounding  rocket  in  the  United  States.  For 
this  reason,  an  imaginary  rocket  with  characteristics  compara¬ 
tively  similar  to  the  Aerobee  is  adopted  to  look  into  the  per¬ 
formance  of  the  sounding  rocket.  Some  characteristics  involved 
in  the  analysis  of  this  hypothetical  rocket  are  listed  below: 

Lift-off  weight,  K  «  1,000  lb 
Diameter,  d  =  13.55  in. 

Fuel  ratio  (ratio  of  weight  of  propellant  to  lifc-off 
weight),  FR  =  0.7 
Specific  impulse,  Isc  =  200  sec 
Drag  coefficient,  CD  -  0,4. 
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With  a  specified  thrust-to-weight  ratio  and  an  initial 
icioe.iLy,  one  can  .integrate  numerically  the  equations  of  motion 
-o  determine  the  maximum  altitude  reached-by  the  rocket.  Thrust- 
ti-’-ae  i gh i  ratios  in  the  range  of  1  through  12  and  initial  veloc- 
i  !■  res  of  0  ,  500  ,  and  1,000  fps  were  selected  in  the  computations, 
i dimputational  results  are  plotted  with  thrust-to-weight  ratios 
ns  abscissas  and  maximum  altitudes  as  ordinates,  as  shown  in 
Fig.  4-1.  The  time  integral  of  the  aerodynamic  drag,  defined  as 


Q  = 


uf 


(Drag)  dt, 


(4-1) 


and  its  relationship  with  the  rocket  performance,  was  also 
s  tudied . 

Figure  4-1  shows  that  the  optimum  thrust-to-weight 
ratio,  (T/W)0pt,  is  located  around  2.7  irrespective  of  initial 
velocities,  and  the  altitude  increases  56  per  cent  due  to  an 
imparted  initial  velocity  of  1,000  fps  at  this  thrust-to-weight 
ratio.  The  slopes  of  the  YInax  versus  T/W  curves  are  much  steeper 
for  values  of  T/W  <  (T/W)  t  in  comparison  with  the  case  when 
T/W  >  (T/W) opt.  The  plots  of  Q,  the  integrated  drag,  indicate 
dut  v  increases  with  T/W  as  expected.  However,  it  is  interest¬ 
ing  to  note  that  the  Q  plots  for  V0  =  0  and  V0  =  500  fps  are 
practically  iuentical  for  T/W  <  5.  This  implies,  '•hat  the  total 
energy  loss  due  to  aerodynamic  drag  is  not  increased  by  impart¬ 
ing  initial  velocities  up  to  500  fps  or  more. 


Best  Available  Copy 


MAXIMUM  ALTITUDE  —  10  FT 


28 


It  is  also  interesting  to  examine  the  percentage 


Computations  have  been  made  fcr  different  thrust- to -weight 
ratios.  Figure  4-2  shows  that  the  percentage  increase  in 
maximum  altitude  increases  nearly  linearly  with  the  initial 
velocity.  It  should  be  indicated  that  the  percentage  increase 
ir.  irsas  is  minimum  at  the  (T/*l)  This  is  no  surprise  because 

Sr 

the  Y_ax  at  the  (T/w) is  much  higher  than  that  at  the  other 


It  is  also  interesting  to  observe  the  percentage  increase 
Ln  payload  with  the  initial  velocity.  It  is  determined  in  two 


steos : 


i)  A  maximum  altitude  is  computed  for  a  rocket  launched 


with  zero  initial  velocity. 

ii)  An  initial  velocity  is  specified.  The  fuel  ratio 
is  adjusted  by  trial  and  error  until  the  maximum  altitude  com¬ 
puted  for  this  rocket  equals  that  computed  in  step  i) . 

The  net  effect  is  an  increase  in  the  payload  delivered 
to  a  specific  altitude.  In  these  computations,  it  is  assumed 
that  the  weight  of  the  rocket  structure  is  10  per  cent  of  the 
weight  of  the  oroceilant,  W  .  It  is  defined  that 
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INCREASE  IN  PAYLOAD 


The  hypothetical  rocket  has 
under  the  assumptions  of  constant  t 
thrust}  for  the  burn  phase  and  a  cc 
throughout  the  vertical  ascent.  I: 
increases  vith  increasing  altitude 
depends  ca  35ach  number-  The  saris 


:StXCatCd 


are. 


srerore. 


underestimated  due  to  the  neglect  of  pressure  ossst  and 


curate  estimate  or  constant  drag  coefficient.  5-r^ever ,  the 
{T/H}  determined  under  these  simplifying  assumptions  deviates 
only  slightly  from;  that  determine  c-y  consissrina  actual  per¬ 
formance,  as  will  be  sees  is  Chapter  V  CSectids  5.3  ass  Fig. 

5-6}  . 

On  the  basis  of  this  investigation,  the  following 
observations  say  also  bz  made: 


i}  The  (T/»)  0p£  seens  not  to  vary  significantly  v£rh 

the  imparted  initi.il  ve  loci  tv.  jn  the  present  case .  Ji/s) , * 

=  2.7  for  y,  =  0  to  V;  =  1,00$  fps. 

ii}  The  decreases  sere  drastically  for  7/K  < 


(T/K)  cpr  ^as  7/*  >  CT/  _ _ 

iii)  The  total  energy  less  sn 


tc  aer- 


increases  with  increasing  T/K,  1st  it  does  not  increase  signi¬ 
ficantly  {=  13%  for  ¥s  =  1,000  fpsl  due  tc  the  inearth  initial 
velocity.  For  the  hypothetical  rocket,  Q,  the  total  drag  inte¬ 
gral,  is  nearly  the  same  for  0  <  ¥s  <  500  fps  ad  T/K  <  5, 


iv)  The  percentage  increase  is  increases  nearly 
linearly  with  the  initial  velocity.  It  is  aiaisss  idas  ®A' 
*  {T/S)c^t.  A  36  per  cent  increase  is  obtained  by  i smarting 
an  initial  velccitv  of  1,000  fss  st 


v)  The  percentage  increase  in  payload  also  increases 
•  V  '  i nearly  with  the  initial  velocity.  It  is  minimum  when 
t  T/W)  0p)_ .  A  42  per  cent  increase  is  attained  by  impart- 
tn  initial  velocity  of  1,000  l.ps  at  (T/W)  n,.,.  . 

v.i  )  The  (T/W)  alters  only  slightly  with  the  fuel 
e.  However,  as  the  T/W  increases  beyond  (T/W)0p^,  the 


•  r responding  Ymax  varies  more  with  T/W  for  higher  fuel  ratio 
'.hat;  for  Lower  fuel  ratio  rockets.  The  closer  the  fuel  ratio 
appt caches  unity,  the  steeper  is  the  slope. 

It  is,  therefore,  considered  feasible  to  estimate 
optimum  T/W  of  real  rockets  under  these  simplifying  assumptions, 
tor  preliminary  design  considerations.  Final  design  parameters 
will  still  have  to  be  based  on  more  detailed  analyses  and  compu¬ 
tations  utilizing  actual  wind  tunnel  test  data.  This  is  attempted 
ar.  a  second  stage  design  analysis  in  the  next  chapter  for  one 
particular  rocket,  the  Areas. 


4.3  Real  Rockets 

A  wide  variety  of  real  rockets15  listed  in  Table  4-1 
were  also  investigated.  The  analyses  were  based  on  the  simpli¬ 
fying  assumptions  of  constant  thrust  (neglecting  pressure  thrust) 
and  a  constant  drag  coefficient.  The  drag  coefficient  was  assumed 
to  be  0.4,  which  is  considered  reasonable,  for  all  rockets  since 
detailed  experimental  data  on  these  were  not  available.  Compu¬ 
tations  were  made  with  zero  initial  velocity  and  thrust-to-weight 

ratios  in  the  range  of  1  to  8 .  Curves  of  Y  versus  T/W  for 

J  max 
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Table  4-1,  Characteristic.* **  of  Real  Rockets. 


Areas  (old)** 


Asp  an 


76.4 


Exos 


Iris 


Meteo 


N; ke-Cajun 


Pol  2 


Skylark 


Veronique  AG-1 


Viking  No.  1 


Viking  No.  II 


Wac  Corporal 


1,500 


5,821 


1,128 


1,500 


1,550 


2 ,560 


4.45 


16.5 


22.88  6.4? 


12.13 


210* 


200* 


0.56 


0.4  * 


3.86 


5.52 


0.4  *  11.1 


230-v  0.8  * 


7.66 


3,040 


9,650 


15,005 


27,376 


2.05  205  0.69  6.47 


*  Estimated  by  the  author,  since  data  not  available 

**  Data  irom  White  Sands  Missile  Range,  New  Mexico. 
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the- so  rockets  are  presented  in  Pig.  4-5. 

'.hi-  foi  owii.u  observations  can  be  made  from  the  plots  of 


1}  For  T/W  •-  2,  tne  performance  of  most  rockets,  in 
r  the  V  attained,  are  degraued  drastically. 

"lac  '  -  1 


or  incline 


dovnwarus  as  the 


x  x ,  ..it  Cv.x ves  vbv>o >£  i 
thrust- to -weight  ratios  increase  bevonc 

ep  t 

iii )  The  i  ; ooe  of  the  curve  for  T/W  >  (T/W)  depends 

opt  17 


(  I  /  «*<  ) 


on  the  1 

n  i  t-i  r  ,'V.nci  ♦“  «  » 

C  Clti  1 . UCiOi.  —  v 

term  (W/d‘)  and 

the  fuel  ratio.  For  the 

rockers 

having 

the  same 

fuel  ratio;  the 

higher  the  W/d2 ,  the  flat 

ter  the 

curve. 

For  ins 

tance,  the  Pol  2 

with  W/d2  =  41.6,  the 

curve  becomes  nearly  flat  as  T/W  increases;  on  the  other  hand, 
the  Meteo  with  W/d2  =  5.02,  the  curve  becomes  inclined  downwards. 
This  results  from  the  fact  that  the  drag  loss  is  much  more  sig¬ 
nificant  than  the  gravity  loss  for  those  rockets  having  low  values 
of  W/d2.  Therefore,  much  care  must  be  taken  in  determining  the 
appropriate  T/W  for  low  W/d2  value  rockets  since  their  perform¬ 
ance  is  altered  significantly  by  this  parameter. 

iv)  Although  some  of  these  rockets  have  low  values  of 
W/d2,  their  curves  still  look  somewhat  fleeter  than  those  for 
some  other  rockets  with  the  same  W/d2.  This  is  because  their 
fuel  ratios,  F_ ,  are  as  low  as  0.5  or  less,  which  cause  an  oppe- 


site  effect 

or.  the  curves . 

However , 

the 

of  W/d 2  and 

F„  is  difficult 

ti 

to  estab 

iish 

pci  X  ciTHG  C  r  X  C  1 

or m,  because  of 

the  com 

lexi 

ar.ee  * 
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I t  shOWi< 


noted  that  all  those  rockets  were  invest! 


Although  the  calcu- 

naxtnu.T.  altitudes  are  not  valid  for  their  actual  performance, 

is  t  ns leered  reasonably  close  to 

opt 

;e  {Sec  Section  5.3  and  Fig.  5--6)  .  It  seems 
that  test,  of  the  rockets  oreviously  investigated  were  well  ae- 
sicr.ed  as  can  be  seen  from  Fig.  4-5.  Looking  at  one  particular 
rocket,  the  Areas,  which  is  used  extensively  for  meteorological 
soundings  by  the  U.  S.  Army,  it  is  apparent  that  the  present 
Areas  is  an  improvement  over  the  older  Areas  design.  By  the 
preliminary  design  analysis,  it  seems,  however,  that  the  (T/W) 
of  this  rocket  is  lower  than  the  present  design  T/W.  Since, 
with  initial  launching  velocity,  the  selection  of  the  optimum 
T/W  becomes  critical  (Fig.  4-1) ,  it  was  decided  to  investigate 
in  more  detail  the  performance  of  the  Areas. 

A  preliminary  parametric  analysis  of  the  performance  of 
Areas  was  made  under  the  assumptions  used  for  other  rockets  in 
this  chapter.  The  results,  regarding  the  influence  of  drag 


_  w  ■ —  —  . 


on  ootimum  T/W;  arid  variation  s  of  Y  ,  V,  and  Y, 

*  luclX  *. 


are  clotted  in  rigs.  4-6  and  4-7,  respectively.  It  is  apparent 
that  variation  in  C_  does  not  alter  (T/W)  significantly, 

^  U  CpC 

pve-  thoucih  Y  would  undoubted lv  be  influenced-  This  esta— 

"  '  *  nieix  "* 

blishes  the  validity  of  assuming  constant  Cri  for  all  rockets  in 
Fia.  4-5  for  determining  optimum  T/W.  Figure  4-7  shows  that 
variation  of  T/W  has  opposite  effect  on  VK  and  YH,  so  that  for 


.  i  .  .  r  m  . *  *  fm  /»•  l  a 

onxv  one  vaxuo  or  -  u/t,)  ,  ~.e 

■*  Cv  =- 


SH3X 


1"  optimized. 
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BURNOUT  ALTITUDE  a  MAXIMUM  ALTI7JDE  - 


i  T/W  «=  (T/W) 


THRUST- TO- WEIGHT  RATIO 

itiorss  or  ournout  velocity  and  burnout  altitude  with 
ratio  for  tbs  Areas,  under  simplifying 
nations  {pressure  thrust  neglected,  c_.  =  0.4,  and 


BURNOUT  VELOCITY  -  SO3  FPS 
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These  plots  for  Areas  were  made  for  most  simplified 
conditions.  A  much  more  accurate  analysis  of  its  realistic 
performance  is  presented  in  the  next  chapter. 


* 


43 


analysis  (Chapter  IV) ,  it  was  round  that  the  present  design  T/W 
of  the  Areas  seems  to  be  higher  than  its  optimum  value.  Since, 
with  an  initial  launching  velocity,  the  selection  of  the  optimum 
T/W  becomes  critical,  it  is  considered  necessary  tc  investigate 
its  performance  in  more  detail. 

In  this  optimization  study,  the  thrust- to-weight  ratio 
is  varied  by  changing  the  propellant  flow  rate  in  the  nozzle. 

This  is  controlled  by  modifying  the  nozzle  size  and  the  exposed 
burning  surface  of  the  grain.  The  geometrical  configuration 
of  the  nozzle,  however,  is  maintained  unchanged.  Other  para¬ 
meters  such  as  specific  impulse,  fuel  ratio,  lift-off  weight, 
dimaeter  and  chamber  pressure  are  also  assumed  to  remain  the 
same. 

The  effects  of  the  pressure  thrust  cue  tc  variations 
of  atmospheric  pressure  with  altitude,  and  variable  drag  coeffi¬ 
cient  as  a  function  of  Mach  number  may  have  been  included  in 
this  study.  Thus,  this  analysis  provides  a  reasonably  good 
analytical  simulation  of  the  estimated  real  performance  of  Areas. 


DRAG  COEFFICIENT 


LEGEND  ‘ 

▲  EXP.  DATA  DURING  BURN  PHASE 
•  EXP.  DATA  DURING  POWER-OFF 


APPROX.  CURVE  DURING 
'  I  BURN  PHASE 


0.3  h 


APPROX.  CURVE  DURING 
|  POWER-OFF 


MACH  NUMBER f  M 


Fig.  5-2. :  Variations  of  drag  coefficient  with  Mach  number  for 

the  present  Areas.. 
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This 


non-cLmensionaiized  sea-level  thrust  history  is  considered 
I  for  other  thrust- tc-weignfc  ratios  and  is  expressed  in  the 


following  simple  equations: 

TSL  '  ?SL.  <0-953  +  1 -03 


for  -  0  <  t/t.  <  0.2, 

—  D 


TSh  TSL  *1*265  ”  0.53  t/tb> 


for  0.2  <  t/tfa  £  0.5, 

Tsl  =  TgL  (1.0543  -  0.1086  t/tfa) 
for  =  0.5  <  t/ 1^  <_  G.933,  and 

TSL  =  *SL  <12*059  '  n*9  . 

for  G.-.933  <  t/t^  1, 


in  which,  TSb  is  the  average  sea-level  thrust  ahd  is  given  by 
the  product  of  propellant  flow  rate  (in  Ib^/sec)  and  specific 
impulse.  .  = 

The  C„.  curves  during  burn  phase  and  during  power-off, 

u 

based  on  the  data  obtained  from  White  Sands,  are  closely  approx¬ 
imated  by  a  number  of  straight  lines  as  shown  in  Fig.  5-2  and 
given  by  simple  equation  as  follows: 


During  burn  phase 


i ) 


CD  ' 

for 

M 

< 

O 

+ 

. 

o 

II 

Q 

U 

i- . 

073  CM  -  0.35 

for 

Zr 

<  < 

1.2, 

Cp  =*  0.321 

-f- 

G.I30S  (3.5  - 

M) 

for 

2 

<  4*: 

3.S, 

CD  =  0.215 

+ 

0.053  {5.5  - 

M) 

for 

-»  * 

V 

<  H  < 

5.5, 

=  0.215 

7 

for 

> 

5.5, 

ii)  During  power-off 


cn 

0.2  +  0.1417  M 

for 

M  < 

1 

CD  = 

0.305  +  0.1075 

(2.5  -  H) 

for 

1.5 

<_  M  <  2.5 

CD  * 

0.25  -f  0.0366 

(4  -  M) 

for 

2.5 

<  'M  ^  4 , 

CD  = 

0-i25 

for 

M  > 

4. 

It  is  obvious  that  the  curve  during  burn  phase  changes 
with  the  propellant  flow  rate.  Since  the  experimental  C0  curve 
is  known  only  for  the  present  Areas  (T/ft  =  4.5) it  is  assumed 
that  the  increment  of  CD  due  to  burning  is  linearly  proportional 
to  the  propellant  flow  rate.  Therefore,  on  the  basis  of  the  CD 
curve  available  for  the  present  Areas,  the  curve  for  other 
thrusts to-weight  ratios  during  burning  may  be  estimated  as 
follows: 


r. 


=  (C  ) 

T/W 


(cD) 


burn  phase 


-  (C.J 


1 


power-off  1 


nnwpr-nff 


so 


'he  approximate  expressions  for  non-diaensionalised 
L  thrust  history  and  curves  as  discussed  above  have 

v 

:  - vq**  coiiuiUtstilo^al  s  tv ,  _ 


h  chang*  ef  thrust- to-veight  ratio  for  a  rocker  ssay 
be  aci :ornpl i shed  by  several  methods*  One  of  the  simple  ways  of 
changing  T/W  for  a  rocket  of  the  sase  weight-  is  by  merely  chang¬ 
ing  the  nozzle  size  and  achieving  a  corresponding  burn  rate. 
Accordingly,  it  is  assisted  here  that  T/K  is  altered  by  a  change 
of  the  nozzle  size.  Further ,  another  reasonable  assumption  that 
is  made  here  is  that  the  area  e  £  the  nozzle  exit  is  linearly  pro¬ 
portional  to  the  propellant  flow  .rate.  Based  on  these  consider¬ 
ations.'  for  the  present  Areas  wit’;  w  =  1  Ih^/sec  and  Ae  =  2.55 
sq.  in.,  if  ?/W  is  changed,  one  can  readily  determine  the 

for  any  T/W  {or  a  corresponding  w) . 


5.3  Optimization  of 7 the  Areas  Rocket 

-  To  maximize  the  obtainable  altitude  reached  by  Areas, 
the  thrust- to-weight  ratio  is  altered  by  changing  its  nozzle 
size  and  propellant  flow  rate.  This,  in  turn,  changes  the  sea- 
level  thrust  history  and  the  curve  during  burn  phase.  On 
the  basis  of  the  information  on  the  present  Areas  {T'/W  =  4.5), 
the  area  of  nozzle  exit,  sea-level  thrust  history  and  CD  curve 
during  burn  phase  are  estimated  for  any  other  T/W  as  described 
is.  Section  5.2.  For  optimization  studies,  these  estimates  will 
be  used  here,  furthermore,  since  the  thrust  of  the  rocket  not 
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only  includes  the  aoaentus  thrust  due  to  rocket  propulsion  but 
also  the  pressure  thrust  due  to  variations  of  altitude,  this 
pressure  thrust  influence  has  also  been  included  for  optimization 
analysis.  The  total  thrust  was  thus  calculates  by  the  previously 
sectioned  relation 
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Numerical  integrations  of  the  equations  of  ration, 
including  the  variations  of  the  thrust  and  the  drag  coefficient, 
were  carried  out  to  solve  the  maximum  altitudes  for  T/N  in  the 
range  of  1  to  6  with  initial  velocities  of  0,  200,  400,  600,  800 
and  1,000  fps  respectively.  Plots  of  Y__„  versus  T/b  for  differ- 

_•  "Goa.  ' 

ent  initial  velocities  are  presented  in  Fig.  5-4. 

The  following  observations  may  be  made  from  Fig.  5-4; 
i)  The  optimum  thrusi-to-weight  ratios  are  practically 
identical  for  different  initial  velocities,  r 


W»>oot  "  2*4‘ 


ii }  For  T/W  <  2,  the  performance  of  Areas,  in  terms  of 

Yroax,  is  degraded  quite  steeply. 

-  iii)  For  T/W  >  (T/W)  . ,  Y  decreases  nearly  linearly 

opt  max 

with  increasing  T/W.  Its  slope  increases  only  slightly  with 
increasing  initial  velocity. 


Iv)  For  each  initial  Velocity,  about  13 


con  c  rncreas* 


in  rs  ostaiEsc  ssea  ?/g  is  siterea  fzo~ 

value  Of  4.5  to  2,4,  which  is  datersinsd  fa-/ 


iSffist  ggg] 


desired  cptlmss  f/ft  for  the  Areas  rochet;. 

v}  The  sain  in  _ due  to  the  is 

ed  by  the  VAKB  system,  is  very  significant 

os r  cent  increase  in  ¥  is  obtained  b"  i 

sax 

velocity  of  1,000  fes  wh=a  T/S  =  fr/sf) 


tne  surnsut  vtiOCita.es  and 


different  initial  velocities  when  T/**  =  CT/S) 

-  #  n/ 

are  uivsn  in  Table  5-1. 


ami  i/ft  =a.5, 


Ttie  integrated  crag,  Q,  throughout  the  f '  * 


■at  as  a  sea- 


sure  of.  the  total  energy  lo'ss  due  to  aerodynamic  drag*-"  Figure 
3-5  shows  that  the  Q  increases  with  the  initial  velocity.  How¬ 
ever,-  the  percentage:  increase  ih  Q  C for  T/&  is  2.4}  dee  to  _x  itial 


velocity  up  to  €00  fps  is  only  6.5  per  cent  or  less. 


t  is  i: 


per  cent  for  initial  velocity  of  1,000  fps.  Thus,  we  conclude 
that  the  increase  in  the  ‘integrated  drag  loss  due  he-  the  impart¬ 
ed  initial  velocity  up  to  600  fps ,  is  almost  negligible. 

It  should  be  noted  that  all  the  computations  for  the 
realistic  performance  estimates  of  Areas  as  given  in  Figs.  5-4 
and  5-6  are  still  based  on  certain  postulations.  For  defining 
the  present  optimisation  problem,  the  gross  weight,  payload, 
structural  weight,  the  propellant  weight  and  its  specific  impulse 
have  all  been  maintained  unchanged  for  the  Areas  rocket.  These 
postulations  were  -felt  necessary  to  sake  the  present  aerodynamic 


performance  optimization  study  tractable.  However, 


>r  a  rear 


design  study,  variations  in  structural  requirements  and  propulsion 
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lesion  (e.g.,  chemical  characters  of  Ihe  propellant,  burning 
-  alterations  and  nozzle  variations)  should  also  be  in- 
tc  was  nor  possible  to  incorporate  these  considerations 
ne  present  Stacy ,  It  is  relt,  however,  that  the  aerodynamic 
limitation  procedure  as  developed  here  is  an  important 
tirst  step  towards  accomplishing  a  final  optimized  design. 

The  performance  analysis  computations  made  for  Fig.  5-4 
were  reasonably  involved  and  complicated.  For  this  reas'"r'. ,  it 
may  be  desirable  to  evolve  a  somewhat  simpler  though  less  accu¬ 
rate  method  of  performance  estimation.  Such  a  method  would  be 
cf  help  in  a  preliminary  design  study.  When  the,  final  design, 
which  includes  the  structural  and  propulsion  considerations,  has 
been  firmed  up,  the  more  accurate  analysis  developed  above  may 
be  used;  ,for  real  performance  estimation.  With,  this  in  mind, 
computations  for  Areas,  have  been  made  with  three  cases  under, 
different  simplifying  assumptions ,  to  see  which  one  is  the  clos¬ 
est  to  the  real  performance  case.  These  three  simplifying  cases. 


are: 


a)  Thrust  constant  throughout  burn  phase  and  CD  varied, 
b.)  Thrust  constant  throughout  burn'  phase  and 
constant  {C^  =  0,43).- 

c)  Thrust  (neglecting  pressure  thrust)  constant  through¬ 
out  burn  phase  and  constant  (CD  -  0.43). 

In  cases  (a)  and  (b)  ,  tl.a  pressure  thrust  is  estimated 
in  advance  and  then  added  to  the  momentum  thrust.  Since  it  was 
found  that  the  atmospheric  pressure  decreases  approximately  lin¬ 
early  -with  time  during  burn  phase,  the  pressure  thrust, 

•aG  0|r,  -  P  )  ,  thus  also  increases  linearly  with  time  appraxmately 
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Once  the  burnout  altitude  is  known,  the  ambient  atmospheric 
pressure  at  this  altitude  is  determined.  And  the  average  pres¬ 
sure  thrust  can  then  be  estimated.  Plots  of  Y  versus  T/W 
for  cases  (a) ,  (b) ,  (c)  along  with  the  actual  performance  as 
presented  in  Pig.  5-4,  are  given  in  Fig.  5-6  for  comparative 
purposes. 

In  case  (a),  (T/w) Gpt  ~  2.6.  In  case  (b) ,  (T/W)0pt  ~  2.8. 
In  case  (c) ,  (T/W)-Qpt  -  2.8.  The  (T/K)  t  value  is  interestingly 
not  altered  much  for  any  of  these  cases.  Thus ,  for  preliminary 
design  considerations,  to- select  the  most  desirable  T/W,  it  is 
considered  feasible  to  achieve  this  by  the  simplest  of  these 
methods,  represented  by  case  (c) ?  i.e.,  when  both  thrust  and 
drag  coefficients^  are  constant  and  the  ^pressure.  thrust  is  also 
neglected.  This  approach  was  used  to  determine  the  optimum  T/W  -- 
for  over  a  dozen  different  rockets  discussed  earlier  (Fig.  4-5)  . 
Results  of  that  analysis  regarding,  optimum  T/W  are  therefore 
considered  valid. 

For  estimation  of  however,  case  (c)  is  not  suitable. 

'  For  Areas,  the  error  in  in  this  case  is  approximately  14 

per  cent.  This  error  is  primarily  due  to  the  neglect  of  the 
pressure  thrust.  It  seems,  therefore,  that  the  contribution  of 
pressure  thrust  to  is  quite  significant  even'  though  the 

pressure  thrust  itself  is  only  a  few  per  cent  of  the  momentum 
thrust. 

Therefore,  for  reasonably  close  estimations  of  both 
(T/W) opt  and  Ymax#  case  (b)  is  considered  suitable.  It  does 
require,  however,  tiiat  the  designer  choose  the  constant  and 
constant  thrust  quite  intelligently.  Estimation  of  the  pressure 
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thrust  contribution  can  be  made  readily  for  rockets  which  have 
a  burnout  altitude  of  60,000  ft  or  more.  .Since,  above  this 
altitude  the  ambient  atmospheric  pressure  is  less  than  7  per 
cent  of  the  atmospheric  pressure  at  sea- level,  the  pressure 
thrust  term  can"  be  estimated  on  the  assumption  that  the  ambient 
atmospheric  pressure  at  burnout  is  zero.  As  discussed  before, 
the  variation  of  pressure  thrust  term,  (PgL  “  pd)  •  i-s  nearly 
linear  with  time-r  So,  a  good  approximation  for  including  the 
pressure  thrust  contribution  for  case  (b)  will  be  to  include 

the  term  A  {P/2)  along  with  the  assumed  constant  momentum 

G  bL  *  ;  "  ‘ 

thrust  for  the  whole  duration  of  the' powered  f  light. 

5.4  Stagnation  Temperature  and  Pressure 

Missile  or- aircraft  flights -at  high  rspeeds  arid  high 
altitudes  introduce  many  new  problems  to  the  designer.  One  of 
the  more  important  problems  that  arises  is  caused7  by  the  skih 
temperatures  (aerodynamic  heating)  that  are  attained  aj:  very- 
high  • velocities .  Adiabatic  wall  temperatures  under  these  con¬ 
ditions  can  exceed  temperature  limitations  of  structural  materi¬ 
als  commonly  used  in  the  manufacture  of  such  missiles.  Although 
the  maximum  velocity  of  the  Areas  is  not  greater  than  4,000  fps, 
it  is  Still  considered  necessary  to  have  the  information  concern 
ing  the  stagnation  temperatures.  -  : 

The  stagnation  temperature,  defined  as  the  total  temper¬ 
ature  of  the  gas  at  rest,  is: 
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’o  =  1‘  +  V-/2C 

a  p 


stagnation  temperature, 

ambient  temperature , 

free-stream  velocity, 

specific  heat  at  constant  pressure. 


It  is  known  that  C  is  not  signif icantlv  affected  by 

temperature  ana  may  be  assumed  to  be  a  constant.  For  air,  it 

is  assumed  that  C  =  C.25. 

p 

Stagnation  temperatures  for  the  Areas  during  optimum 
performance  as  computed  from  the  above  relation  are  shown  in 
Fig.  5.-7  and  the  maximum  values  for  different  initial  velocities 
are  1.  '•  in.  Table  5-2 .  . 

It  should  be  indicated  that  the  maximum  stagnation 
temperature  is  reduced  when  the  T/W  is  altered  from  its  present 
design  value  (T/W  =4.5}  to  that  of  the  proposed  optimum  case 
<T/W  =  2.4}.-  In  fact,  for  initial  velocities  of  nearly  500  fps, 
the  maximum  stagnation  temperature  for  the  optimized  case  Is 
below  the  maximum  stagnation  temperature  for  the  present  Areas . 
Even  with  an  initial  velocity  of  1,1)00  fps,  this  temperature  for 
the  optimum  case  is  only.  1,680  WR  Which  can  be  readily  withstood 
by  the  present  nosecone  materials.  Similarly,  the  external  skin 
temueranures  are  expected  to  ba  within  reasonable  design  limits. 

It  can  thus  be  concluded  that  the  aerodynamic  heating  of 
Areas  due  to  Importation  of  initial  velocity  is  no  problem  for 

6K1.S  ^2  iiC*  * 
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The  stagnation  pressure,  defined  as  the  sum  of  the 
ambient  and  the  dynamic  pressures,  is  given  by 

Pc  =  +  dV2/2 

a 

where 

Po  =  stagnation  pre'ssure, 

P_  =  ambient  pressure,  . 

p  =  ambient  density,. 

V  =  free  stream  velocity. 

Variations  in  stagnation  pressure  during  powered  flight  for  the 
present  Areas  and  for  the  proposed  optimised  Areas  are  presented 
in  Fig.  5-8.  It  can  be  seen  that  the  maximum  stagnation  pressure 
is  reduced'  from  25  psi  to;  16  psi  when  the  T/K  is  altered  from  its 
present  design  value  to  the  proposed  optimum  value.  Even  with  an 
imparted  initial  velocity  of  1,000  fps,  the  maximum  stagnation 
pressure  (23  psi)  for  the  optimum  case  is  still  lower  than  for 
the  present  Areas.  This  is  well  within  the  existing  design  limits 
and  thus  represents  no  problems. 


STAGNATION-  PRESSURE 


CHAPTER  VI 


CONCLUSIONS  AND  SUGGESTIONS 
FOR  FUTURE  WORK 


Significant  gains  in  performance  c-f  single-stage  rockets 
are  considered,  possible  by  utilizing  the  VAHB  system.  The  gain 
in  maximum  altitude  of  an  optimized  Areas  when  launched  with  an 
initial  velocity  of  1,000  fps  is  expected  to  be  93  per  cent  above 
the  estimated  maximum  altitude  of  the  present  design  Areas  with¬ 
out  bobs ter.  comparable  gains  for  other  single-stage  rockets 
are  also  expected.  It  seemsr  that  single-stage  rockets  designed 
in  an  optimum  fashion  for  conventional  launch  will  continue  to 
perform  near  optimum  even  with  the  VAMB  system.  This  was  found- 
r  to  be  the  case-  for  Areas.  For  other  rockets,  however,  a  similar 
optimization  study  should  be  performed.  It  is  felt  that  beyond 
the  aerodynamic  optimization  as  presented  in  this  thesis,  an 
integrated  optimization  study,  for  the  individual  rocket,  includ¬ 
ing  the  structural  and  propulsion  design  variations,  should  be 
performed  to  evolve  the  final  optimum  design. 

It  is  feasible,  for  preliminary  optimization  study,  to 
assume  that  both  thrust  and  coefficient  of  drag  are  constant  and 
the  pressure  thrust- is  neglige -le.  The  estimation  of  optimum 
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PROGRAM  1 

l PAST  SQUARES  CURVE  FIT  FOR  VARIATION  OF  ATMOSPHERIC  DENSITY 
rill k. ALTITUDE 


D I MENS  ION  A ( 3 , 3  5 , R ( 3  5 , Y ( 20 ) , RHO  ( 20 ) ,RHGC ( 20 ) , Z ( 20  > 

RFADH»I)  K-r N 

1  FORMAT! 2  I? 5 

READ!  1,2)  (YU  )  »RHO(  I  )  , I=1»K)  ' 

2  FORMAT! 2F10.0) 

Ai'l,2)  =  0.  ' 

A  (  i  «  3  )  =  0 . 

A  {  2 ,  3  )  =  0 » 

A (3,35=0.  ~ 

8  ( 1 )  =  0  » 

B ! 2 5 =0 . 

8(35-0. 

00  io. 1=1,12 

A ( 1, 2  3  - A { 1 , 2  >  +Y (  I  )  ' 

A(1,3)  =  A<  W3  5+YU  5**2 

A(  2 , 3)=A.(  2 ,3  )+Y(  I  )#*3  -  -  - 

-■  A  {3:,  3  5=  A  ( 3  » 3:)'-l-Y-:{ !;) **4  ' 

8(1)  =  8  f  1  5  +  RHO.M  j  -  -  v 

:  3  (  2>=8(  2-)  +  RH0  m*YUT) 

8!  3')— 8(  3  )+RHO  (ST  5 *Y  { J  5**2  .  . 

TO  CONTINUE  “  ‘  _  _  . 

r  AM,  15  =  12. 

A ( 2 , T) =A ( i , 2  5 
'  A( 2,2)=A( 1,3} 

A  C  3, 1 ) =A ( 1 ,3  5 
A( 3, 2 )=A ( 2 , 3 ) 

CALL  S  IM0(A,8,N,N,KS) 
jF(KS)  SOM  1 ,  30 
11  WRITE ( 3  » 1 2  5  B(1  ),B(2 5,B(3i 

>2  FORMAT?  1H1 ,2X,3HB1  =  ,E16.-R,5X,3H62=>E16*8,5X,.3HB3  =  ,E16.8  ) 

WRITE  (3,-13  5 

13  FORMAT  (  10X  ,  8HALT I TUDE  , 10X  * 1  OH  ACTUAL  RHOtlOX, ],  5H  APPROX  I  MATE  RHO) 
DO  20  I  =  U12 

RHOC( 1 1 - ?  1 1 >+B( 2  5*Y  M )+B(3)*Y(I 5**2 

20  WRITE! 3,21 )  v  < I ) » RHO ( I ) , RHOC ( I ) 

21  FORMAT { 5X , F i 0 . 1 , 8X  y  E16 . 8 , 6X , E 16 . 8  5 
SUMY=0. 

Sum 2=0. 

SUMY /=0. 

SuMY2=0. 

DO  15  :=12,K 

Z( ! )=AL0G10(RH0(I ) ) 


SUMY=SUMY+YU  ) 

SUMZ=SUMZ+ZU  J 
SUHYZ-SUMYZ+YU  )*Z(  I  ) 

15  SUMY2=SUMY2+Y( 11**2 
D=5L0AT(9)*SUMY2“S0MY#*2 
DA=S,JM7.  ^Sl/MY2“'SUMYJ*tSUMYZ 
DC=FLOAT ( 9 )*SUMYZ-SUMY*SUM2 
P=10.**(DA/D) 

C=OC/D/AtOGiO(?  t?18 ) 

DO  25  1^12  *K 

RHOC?  U  =  P*EXP(C*Y(I  } ) 

25  WRITE(3,26>  YU  )  fRHO  (  I)  tRHOC  (I  > 

26  FORMAT? 8X.F10.1 ,8X»E16.8f 6X.E 16.8 ) 
WRITE { 3»2? )  P,C 

30  ST0DAT1 1X,?Hps*6i6-8*5Xt2HC=f E16.8) 
END 

iDATA 

20  3 

0. 

5U00, 

loobo. 

15000. 

20000. 

25000. 

30000. 

35000. 

40000. 

45000. 

50000. 

55000. 

60000. 

65000. 

70000. 

75000. 

80000. 

8506g. 

90000. 

100000. 


0.07647 
0.06590 
0.05648 
0.04814 
6.04077 
0.03431 
.02866 
0.02375 
0.01890 
0.01487 
0.01171 
0.00922 
0*00726 
0.00572 
0.00448 
0.00351 
0.00276 
0.00217 
'  0,00171 
0.00107 
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PROGRAM  2 

NUMERICAL  SOLUTION  OF  THE  EQUATION  OF  MOTION  FOR  4  VERTICAL 


DIMENSION  V { 2  I  t DERY *2 )  ? YN 12 Lf YNN { 2 ) » YNNN ( 2 i 

COMMON  AR&A »CD*RHG1 ♦ ALPHA*  WFR,WF,BRATE*87,G, A,6, C  »ANUM»  RHO 

l  «AMACH>TlNF»AltBItC! »P1 *ALPftAi,THAV*P 

:  RH0i=0. 1 29184  - 

5LPHA=-0.480263 E-4 

A=0.757S78F-1  - 

'  P-=’~0.202365E-S 

C=0  .  L50855E.-1O  :  - 

P1=0.2&5263E4 _ _  _ 

ALPHA  l=-r9.4?3iC9£-4 

A1=0-20B395E4  “  _  “  , 

81— -Q«666?94£~l  S 

Xi  =  0.59616SE-6  ^  ' 

.  0=32.2 

:  DT=0.5  •  -I-  -  .  •  '  -  .  . 

00.  REAOC-*  »9)  WFR»0»SP-IM».WF  - 

9  FORMATS  4Fl<T.O )  -  .  ■-  -  v  '  '  .  ' 

veXH=SPIM4G  r  -  .  _ 

■  AREA=d.7;654_-D?D/l»4.  .  -  ' 

-  FRAT'IOil/FAWFR  .  ..  ' '  v  >  r  ...  7 

VZERO=G.  .  '  . 

do  40  J~T-»6  ‘ 

5  -TW-T.  '  :  =  "  - 

00  30  i!  =  lfTil  .  : 

THAV=tw*WFR  -  ■= 

:  8RAT£=Tw*£4WF.ft/VEXH  '  . 

-:bt-wf/8rate 

Y C I ) -V  ZERO 

Y(  2  J=0*;  :  :  .  r  -  1 

T=Oe  - 

‘  HR  LT  E  { 3  *  I J  - 

1  FQRMATi.///  j 2 X ,  lOHTHRuST/WT *  , 2X ,6HFRAT ID »  2  X  *  9H73TA  L  W  T.  ,2  X, 4HFUEL  , 
l8Xt  5H8R'ATE-,4Xt5HBTlME»5Xt7MEXtiVEL*-*I4X*--4HAREA.5X.*5HVZER0> 

WRITE  C  3,2)  T W* F R A T:I0 *  WER*  , BR ATE  » BTT ,  y  £?C'H ,  AREA  t  VZ  ERO 

2  ’FORMAT!  lH0»7F10.3iiq£'t2FiO.3> 

HR  I  TEC  3,3.)  -  - 

3  FOR MA T ( 1 HO , 5X , 1 OHT I HE  f  S  EC  S ) , 5  X.16HVE LOC  I  TV { F  T /SEC ) , 5X , 

1 1 2HALT I  TUBE  (  FT  ) ,  3X  , iOHDRAG  GOEF.  :12X  ,5H0RAGF) 

ORAGr^O.  •  r 

4  CALL  FCT { T ,Y»DERY )  -  : 

fiRAGF- ORAGF'S-ANUM^OT/G  - 

R1=0ERY ( 1 )40T 
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S1=0EP.Y!2)*DT 

TN=T+0T/2« 

YN !  1 )  =  Y  ( 1 )+ R I  /  2  . 

YN( 2)=Y( 2 )  +  Sl/2  * 

IF{YN{ 1) .LT.0.«AND.Y<2).GT.0.)  GO  TO  21 
CALL  FCT ! TNfYNtOERY ) 

R2=OFRY'(  1  )*DT 
S2=DERY (2 )*DT 
YNNt 1 j  =Y{ 1 J+R2/2  « 

YNN!2)=Y!2)tS2/2. 

IFCYNNt 1).LT*O..ANO*Y{2)*GT.O. )  GO  TO  21 
CALL  FCT { TN *¥NN  *D"ERY ) 

R3=DERYU)*DT 
S3=06RY(2»*0T 
_  TNN=T+OT 

YNNNU)=Ym+R3  - 

YNNN  (2 )  =Y  ( 2 1-+S3  -  : 

•fF!  YNNN  M  i  .LT  .O..AMO.Y  12  >  .GT.O.  )  GO  TO  21 
CALL  FCTITNNtYNNN,OERY) 

R4-DERY  C 1 )*0J 
S4=DERY!2)*DT 

DYi=! R1+2.*R2+2.*R3+R4  1/6. 

DY2=<Sl+Z.*$2+2.*S3*S4>/6. 

Y  M)  =  Y  {  1  ) +0  Yl  ' 

YI2>.=Y<2>+0Y2 

T=TNN 

..  I  Ft  A8SI  '7-BT:)'-i  *  1  80*80*81-  -- 

80  ACC  §=6Yl/D7/G  . 

WRITE!  3*111  TNN ♦ Y ( 1 1, Y ! 2 ) «  C  0  * DR  A&F ♦  ACCE 

11  FORMAT!  lH0t4X *F8*2,4X,2:f  5 X *E14.7 )  ,  8X  ,F  10.4*8X *  E14.7  ,F  10*3 

81  IFTT.GT.BT-.AND.YU  >. IE. 0.  )  GO  TO  21 
'  .GO  TO  4 

21  WRITE!  3*221  TVYU  1,Y!2I*DRAGF 

22  FORHATC 1H0*4HT1HE»FT,2*7HSEC0NDS*5X,8hVEL0CITY,E14,7, 
16HFT/SEC.5X*8HALTITUD6*E14.7«2HFTt5X*6HDRAGF=,F14.7J 

-  TW-TW+0.5  - 

50  CONTINUE  - 
VZER0=VZER0+200. 

40  CONTINUE 
101  STOP 
END 

SUBROUTINE  FCT!T*Y,OERYi 

01 MENS ION  Y 1 2 ) ,DERY ( 2 ) * YN ! 2  J , YNN ! 2 ) *  YNNN f 2  i 
COMMON  AREA tCOfRHOl* ALPHA*  WFR*WF*BRATE,BT,G,A,&,C ,ANUH 

1  *AMACH,TlNF*Al*Bl*£l*Pl*ALPHAltJHAV*P 

IF! Y! 2 1 .LE. 50000* )  RH0=A+B*Y!2 )+C*Y< 2>**2 
IF! Y(2J .GT. 50000. I  RH0=RH01*EXP! ALPKA*Y!2) } 

IFCYT2 ). IE. 50000. >  P=A1+61*Y!2 >+Cl*Y! 2)**2  c 

IF! Y! 2). GT. 50000.1  P=P1*EXP< ALPHAI*Y!2  ) ) 

IF! Y!2).GE. 120000.  .AND.  Y ! 2 ). LE. 164000. }  TI NF=45d.+ !Y ! 2 ) 


RHO 


120000.) 


■ 


*18./  4400. 

IFIYC2).G1 -1-64C00.  .AND.  YI2  ).L£. 200000.  }  TJNF=630. 
i  *  I  v  i  ?  $  »GT  *200000.  .  A  NO  .  ¥  i  2  ) .  LE.  2  50000.  )  TiNF=450*+  i  250000*  -Y{  2) ) 

l  UUu. 

IF(V{2i.LT. 120000.  .OR.  Y(2)*GT»250009. )  7INF=p/RHQ/ 53.3 
K  -  49.02*  SORT  O  INF- ) 

AMACH=v{; ) /SONIC 
IPiT.OT.f’Tj  GO  TO  80 
I P i AMA CH .16.0.9)  00=0.3 

I  6i AKACH . GT *0 *9 .AND  , AHACH .LE  • 1.2)  CD -G. 2+1 .073= {  AH*CH-G.9) 

I  =  <  aMACH.GT. 1 .2 -AND. AMACH.L E . 3. 5 )  CD=0.  32 1+, 1309* l 3. 5-AHACHl 
i?[  AMA-fH.OT  ,3.5.AN0.AKACH.LE.5i5)  00=0.215-*-.  053*  C5.5-AHACH) 

IP? AMACH,GT.5„5)  C 0=  . 2 1 5 
001=00 

1 F I AHACh.L E. 1*5)  C0=.2+.i41 7=AHACH 

I F I  A  MAC  H  »  GT . 1 . 5 . AND . AMAC  H .L  E . 2 . 5 )  SO = *  305 +  .  1 07  5*  {  2 . 5-  AHACH ) 

I F  l  AM  ACH  .  G  T  *2 . 5  .  AMO  .  A  MACH.  L  E  .4  .  )  €  0= . 25  +  *  0366*  C  4  .  -AM  ACB ? 
t F { AMACH.GT .4 . )  CD= . 25 
IFIT..GT.BT)  GO  TO  60  . 

CD2=CD 

C  D=C  02  + CCDJ-C02 J*RRATE 
ANUM=0. 5*CD*AREA*RHQ*Y ( 1 ) =*2 
IP*  7—  RT)  61 .81 1 62 
T7=T/f?T 

IF{TT.L£.0-.2)  TSL=THAV*{0. 953+1. 03*7T) 
i F { 77 ,G7 . 0 . 2 . AND . f T . LE . 0. 5  >  TSL=THAV*  C 1 . 265-0.53*77 ) - 
IrCTT. GT. 0.5. AND. TT.LE. 0.933)  TSL=7HAV*I 1.0543-0. 108 6*77} 

I F  C  77  .GT  *0 . 933  wAND .  TT  .  L  C* .  1 . 1  TSL=THAV*  {12.  059-1 1 .9* TT  I 
770TAL-TSL-s-2.55*{2116.-P)/146.*8RATE 
fw=T i GTAL/H-R 

9  £  N  0  m  =  w  F  S-  8  R  A  t  g  *  7  '  - 

DRA6=AfiUM/D5NCH 

06  RY  { 1 )  =7R*G/  1.1 SRATE*7/WFS  )-DkAG-G 
OERVT 25=Y( 1 ) 

[ P ( 05RY Cl? .LE.O.O. ANO.Y (2 }*LE. 0. )  GO  70  71 

RETURN  _ 

OrkY  { 1  }=0. 

Rr  T IJ  F-  fi 
t-BN  QM*WFR— WP 

or  AG= a hum/ den m 

DERYf 1 )=- DRAG-6 
D£kYC2)=Y(1) 

I F{ OERY ( 1 ) .LE.0.0.AND.Y12  J.LE.O. )  GO  70  71 

RETURN 


i*  ,45 


315. 


(74) 


75 


10.  Leitmann,  G.,  "Thrust  Programming  for  High  Altitude 

Rockets”#  Aeronautical  Engineering  Review#  Vol.  16, 

No.  6,  pp.  63-66  (1957) . 

11.  Tsien.  K.  S.#  and  Evans#  R.  C..  "Optimum  Thrust  Programming 

for  a  Sounding  Rocket" ,  Jet  Propulsion,  Vol;  21,  No.  5# 
pp.  99-107  (1951). 

12.  Goddard,  R.  H.#  "A  Method  of  Reaching  Extreme  Altitude" , 

Smithsonian  Misc.  Coil.#  Vol.  71,  No.  2  (1919)#  and 
American  Rocket  Society  (3946). 

13.  Scarborough,  J.  3.#  "Numerical  Mathematical  Analysis”, 

5th  ed. ,  The  Johns  Hopkins  Press,  1962. 

14.  Shapiro#  A.  H. #  “The  Dynamics  and  Thermodynamics  of 

Compressible  Fluid  Flow1*#  Vol.  1,  The  Ronald  Press  Co., 
1953.  ' 

15.  Ordway,  F.  I.,  Ill#  and  Wake ford ,  R.  C. ,  " International 

-  Missile  and  Spacecraft  Guide”#  McGraw-Hill  Book  Company# 
Inc.,  19<j0^ 


?  ©#?*  ~«N  a  TiK£  AC  Tiv*  T  v  iCorpijrsiz  authot)  1 

I  i 

1 

l#.R£PO«T  SECURITY  CLAfStFIC  ATtOW 

Unclassified' 

U.  S .  Army  Research  Ufflee-Burham 

2b.  CROUP 

HA 

4.  OCJCmCTI  VE  NOTES  f  JVfi#  «f  rfpctl  tful  lactUMlw  dtt**) 

Technical  Report  No.  9 


S-  autho«{5»  (Fittt  ci<£{£*  tniti*lr  Ust  n+m+j 

Sudhir  Kumar 
Stanley  Chin-nan  Shi eh 


•-  REPORT  DATE 

August  19o9 

M.  CONTRACT  OH  CRANT  HO. 

*.  RROJCCT  NO- 

S-7-S 

f. 

- 

a. 

TOTAL  NO.  or  RA«CX  [7*.  mo.  or  rcrs 


KM*  npoti) 


IV.  VitriUBUnON  3TAI&ML*I7 

This,  document  has  been  approved:  for  public  release  and  sale;  its 
distribution,  is  unlimited .  - 


It.  SUSSLIh£nT4Ht  HOTO  «.  SRONSOmHC  MILITARY  ACTIVITY  ~  '  ~ 

If.  5.  Army  Research  Office-Durhas 
Box  CH~  Duke  Station  - 
Durham,  North  Carolina  27706 


»*.  ABSTRACT  ~  ‘  -  -  -  -  ■ 

Studies  on  optimization  (altitude  maximization)  of  a  single-stage  rocket 
performance  with  and  without  the  use  of  the  Vacuum-Air  Missile  Boost  (VAMB) 
system  are  presented.  Performance  estimation  and  optimization  of  a  hypothet¬ 
ical  1,000  lb  rocket  with  ar.d  without  the  VMS  system  have  been  done.  Similar 
studies  have  also  been  performed  for  eighteen  resets  of  various  countries' 
without  the  YAM3  system.  It  was-  found  that  the  optimum  tfcrust-to-weight  ratio 
(T/w)  for  conventional  rocket  can  be  estimated  roughly  for  preliminary  design 
of  rockets  under  simplifying  assumptions .  Detailed  performance  and  optimiz¬ 
ation  studies  have  been  presented  on  the  Areas  rocket.  It  was  found  that  the 
optimum  T/W  for  Areas  with  and  without  VAIS  system  is  approximately  2.k,  which 
is  almost  half  of  the  present  design  value  (4.5).  The  gain  in  maximum  altitude 
-  of  the  optimized  Areas  launched  with  an  initial  velocity  of  1,000  fps,  over 
that  of  the  present  design  Areas  without  booster,,  was  estimated  to  be  nearly 
93  per  cent.  - 
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